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Abstract. The behavior of the critical function for the breakdown of the homotopically non-trivial 
invariant (KAM) curves for the standard map, as the rotation number tends to a rational number, is in- 
vestigated using a version of Greene's residue criterion. The results are compared to the analogous ones 
for the radius of convergence of the Lindstedt series, in which case rigorous theorems have been proved. 
The conjectured interpolation of the critical function in terms of the Bryuno function is discussed. 



A long-standing problem in the study of quasi-integrable Hamiltonian systems is the characteriza- 
tion of the threshold for the break-down of KAM invariant surfaces in terms of the arithmetic proper- 
ties of the frequencies vectors. In this context, we consider a simple, yet paradigmatic, discrete-time 
model, the so called standard map, introduced originally in ifTTTI PT~S1 - The standard map is the dynam- 
ical system defined by the iteration of the map 



on the plane R 2 given by the same formula as ( 11.11 ) with (£, 77) replacing (x,y). For some background 
information, we refer the reader to the enormous literature on the topic, and in particular to l2"31 for a 
review. 

Despite its apparent simplicity, there are only a few properties of the standard map which can be 
considered really well understood to full extent, especially from an analytical point of view. For 
instance the existence of KAM invariant curves, for values of the parameter e small enough and 
Diophantine rotation numbers, has been proved a long time ago, but only recently the dependence of 
the radius of convergence on the rotation number has been obtained [ 14 5 1 as an interpolation formula 
in terms of the Bryuno function (see below). Also for the studying of the separatrix splitting, only 
recently the original program by Lazutkin [20| has been completely achieved in a rigorous way fTTll . 

In particular no rigorous analysis has been implemented for detecting the critical value of e at which 
the KAM invariant curve breaks down, and only numerical results and heuristic theories exist on that 
subject; see I23II25I IT1. 



1. Introduction 




(1.1) 



Here (x, )i)eTxl; but of course the map T £ could be lifted to a map 
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In dl.lt we can eliminate the y variable by writing the dynamics "in Lagrangian form" as a second 
order recursion: 

^Xfi H~ x-fi—\ — £ sin jCji , (1.2) 

for all n € Z. 

For s = 0, the circles y = (const.) are invariant curves on which the dynamics is given by rotation 
with angular velocity co = y/2n; we call co the rotation number. Without generality loss we can 
choose co e (0, 1) as the invariant curves of the standard map are invariant under translation of 2n in 
the j-direction. 

As the perturbation is turned on, we face the classical KAM problem of determining which invariant 
curves survive and up to which size of the perturbative parameter s. Such invariant curves are given 
parametrically by the equation 

\x = a + u(a, £, co) , 

\y = 2nco + u(a, e, co) — u(a — 2nco, e, co) , 

where in the a variable the dynamics on the curve C £ , w is given by rotations a n+ \ = a„ + 2nco (which 
solve (II .2I > for e = 0). The function u(a,s,co) is called the conjugating function or linearization, and 
satisfies the functional equation 

[D^.u ) (a, e, co) - u(a + 2nco, e, co) - 2u(a, s, co) + u(a - 2nco, e, co) 
v ; (1.3) 
= e sin(a + u(a, e, co)) , 

whose solutions are formally unique if we impose that u(a, e, co) has zero average in the a variable. 
Therefore the study of the invariant curves C £ , w and of their smoothness properties may be reduced to 
the study of the existence and smoothness of the solutions of the functional equation (11.3b . 

The solutions of (11.31 1 can be studied perturbatively by formally expanding u(a, e, co) in Taylor 
series in e and in Fourier series in a; the resulting series is what is traditionally called the Lindstedt 
series: 

oo oo 

u(a, s, co) = V e k u [k \a, co) - Y e k ^ ^ ■ (1.4) 

k=\ k=\ \v\<k 

To characterize the breakdown of an invariant curve C £iW we introduce the radius of convergence of 
the Lindstedt series 




(1.5) 



the lower (analytic) critical function 

s c (co) = sup{e' > : Ve" < e' C E »^ exists and is analytic}, (1.6) 

and the upper (analytic) critical function 

s c (co) = infje' > : V e" > s' C £ " >a) does not exists as an analytic curve}. (1.7) 

In general one could define analogous functions for negative values of s; for the standard map they 
would be anyhow identical (by symmetry properties). 



SCALING OF THE CRITICAL FUNCTION 



3 



Clearly p(a>) < e c (to) (in the early papers on the subject some confusion was often made between p 
and £ c ). It is instead believed thai, for the standard map, e c (to) = s c (to), so we can speak generically of 
one critical function s c (a>) without further qualification. Note that for similar maps with more general 
perturbations numerical results |3lj suggest that the two critical functions may be indeed different. 
Note also that one could define breakdown thresholds with the analyticity condition in dl.6t . (I1.7t 
replaced by a weaker one (such as C°° or C k ); again those thresholds could, in principle, be different 
from the analytical one, though for the standard map it is believed that no such difference exists, so 
that the analytic category is the right one to investigate the breakdown phenomenon. 

The radius of convergence of the series ( 11.41 ) is zero - and no KAM invariant curve exists - when 
to is rational. When to satisfies an irrationality condition known as the Bryuno condition (see below), 
instead, it can be proved that p(a>) > - so that analytic invariant curves exist for s small - and 
even precise upper and lower bounds on the dependence of p{to) on to can be given, up to a bounded 
function of to fT4l l5l. More precisely for any rotation number to one can define the Bryuno function 
B(to), as the solution of the functional equation 1521 

B(to) = - log to + toB{to~ l ) for to e (0,1) and irrational, 

(1.8) 

B(to + 1) = B(to). 

By an easy fixed-point argument it can be proved that a solution to ( 11 ,8b exists and is unique in L P (T) 
for each p > 1. 

We shall call Bryuno number a number to satisfying the Bryuno condition B(to) < oo. Then for any 
Bryuno number to one has 

\logp{to) + 2B(to)\<C Q , (1.9) 
for a universal constant Co, that is for a constant Co independent of to (see 1141 Bl for a proof); in 
particular this implies that an invariant curve with rotation number to exists if and only if to satisfies 
the Bryuno condition. Equation ( 11.91 ) and similar formulas are referred to as "Bryuno's interpolation 
formulas". 

The claim is often made that a formula analogous to (11.91) should hold for s c (to): for any to satisfying 
the Bryuno condition one should have 

|loge c (w)+j5fl(&>)|<Ci, (1.10) 

for a universal constant C\, with an exponent B < 2; it is conjectured that B - 1 (see e.g. l28l[m K 

Equation d 1 . 10l > implies a scaling law for the critical function s c {to) as to — » p/q on suitable se- 
quences of Bryuno numbers. In fact, given d 1 . 1 Ob . there are sequences of Bryuno (even Diophantine) 
rotation numbers chosen in such a way that along them the critical function tends to zero in any arbi- 
trarily fast way. For example, we can consider the two sequences of Diophantine (even noble) rotation 
numbers 

1 1 

to k = - , to k = , (1.11) 
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where y = (V5 - l)/2 - [1°°] is the golden mean; then (ll.lOt . with /3 = 1, would imply that 
s c (wfc) = 0(\/k) while £ c (a>fc) - 0(e~ k /k), that is much faster (see Q, p. 625-626). It is therefore 
essential to have a good control over the arithmetic properties of the rotation numbers one considers 
when speaking of scaling properties of the critical function £ c (w). 

The conjecture of Bryuno's interpolation was actually made for the critical function e c (ca) more 
than 10 years ago in [28 1; in that paper, Bryuno's interpolation is stated formally for the radius of 
convergence, but the numerical calculations, with Greene's method, compute the critical function. 
The main motivation behind ESI was the comparison with the work of Yoccoz in P2"l . together with 
the claims of universality coming from the renormalization description of the critical invariant circle. 
In ill?! |3 (see also 0) Bryuno's interpolation for the radius of convergence was indeed proved; 
the mechanism of proof in Q, based on the multiscale decomposition of the propagators in the tree 
expansion, naturally generates an estimate of p(a>) in terms of the Bryuno function for the semi- 
standard and standard maps. On the other hand, there is no compelling a priori heuristic reason for 
the critical function e c (a>) for the standard map to satisfy an interpolation formula in terms of the 
same arithmetical function as the radius of convergence p(a>); put it in another way, it is by no means 
obvious that | log p(a>) - (2//J) log £ c ((o)\ should be bounded. 

From this point of view, it would be interesting to consider generalized standard maps, i.e. maps 
where the nonlinear term in (11.1b is an arbitrary analytic, periodic function of x (see 0|6]|). In these 
cases a Bryuno's interpolation formula for the radius of convergence of their Lindstedt series is not 
known. 

The method used in EJ cannot be pushed so far to get reasonable numerical data on the critical 
function, for some rather obvious reasons; in fact, the method we used there (Pade approximants) 
attempts at modeling the whole natural boundary, giving particular weight at those regions of the 
boundary where the singularity is "strongest": that is, to those regions closer to the origin (the first 
order or dominant singularities as defined in EJ), which determine p(cS); so that part of the natural 
boundary near the real s axis, which determines e c (oj), is represented, as a> is closer and closer to 
a rational value, as a few scattered points from which no reliable information can be extracted: this 
happens already for rotation numbers as little close to a rational value as, for instance, 1/(50 + y) is to 
0, that is still quite far from the rational value. One clearly needs a method in which all the computing 
power is dedicated to the calculation of the quantity one is interested in, that is e c (oS). 

To this aim, two methods have been used previously in the literature: Greene's method (also known 
as residue criterion; more about it in the next section), used in [28], and the frequency map analysis 
fl9i used in El- As we also use Greene's method, we shall postpone a more thorough analysis to 
the next section, and go on to a discussion of the results of fill . 

In fTTIl the following functions are defined: 




(1.12) 



SCALING OF THE CRITICAL FUNCTION 



5 



and 

f P i) 

o) p / q (e) = sup |a» < — : C e ,u) exists and is C j . (1-13) 

The meaning of those functions is that, for the given value of e, no (C 1 ) invariant curves exist 
with rotation numbers between oT^(e) and a)^ q (e). The frequency map analysis method computes 
^plq( E ) - a, p/ t y( £ ) _ w p/ ? ( e ) f° r se l ec ted values of e; A p / q (e) should tend to with e and in this way 
a lower bound on fi should be obtained (see below). Note that s is fixed, and correspondingly some 
rotation numbers are computed numerically, therefore losing any strict control over their arithmetical 
properties. 

We remark that the regularity properties of the functions ar^^e) are quite hard to understand, and 
in particular their relation with the critical function is far from obvious. In fact, while it is certainly 
true that 

s c (tb) <s Vwe {o)~ lq {s), (Op /g (e)), (1. 14) 

the formulas at p. 2037 and p. 2052 of ATI , that is edco^^e)) = e, cannot be claimed in full rigor 
since an invariant curve with rotation number very close to p/q can be broken by the effect of another 
resonance p' /q' ~ p/q, but distinct, so that we can at most claim that 

E c {0f pjq {E))<S. (1.15) 

This implies that the law |a»j^ (e) -p/q\ ~ s q , numerically determined in fffl . provides for the critical 
exponent an estimate from below of the actual value fi, which in principle could be higher (if it does 
exist at all). Equality in d 1 . 15l > can be safely assumed at best for e such that the corresponding value 
ojpj^s) belongs to a special class of rotation numbers tending to p/q (in some sense the "best ones", 
that is the ones whose partial quotients grow as slow as possible), which are indeed the ones considered 
in fTT1 and in the present paper (and which are the only really accessible to a numerical investigation). 
Note also that to saturate d!.15t one should assume other qualitative features (like monotonicity) on 
the functions oj p ^ (e), which are far from being proved. However for the case a> — > only this is 
enough, since estimates in [30] imply an upper bound on the critical exponent, which closes the gap. 

The numerical lower bounds found in fTTI for fi are consistent with fi = I with errors of orders 4% 
for d close to 0/1, 10% for oj close to 1/2, 5% for to close to 1/3, 10% for oj close to 1/4, 8% for oj 
close to 1/5 and 10% for a> close to 2/5. 

Establishing a condition like (ll.lOt is out of reach from the numerical point of view if one wants to 
take into account arbitrary sequences of Bryuno numbers. In fact for the frequency map analysis this 
is a limitation intrinsic to the method itself, since it automatically sort of chooses the best sequence 
of Diophantine numbers tending to any given rational value. For any other method, like Greene's 
residue criterion, to investigate Bryuno non-Diophantine numbers would require computer resources 
far beyond current availability, while computing the critical function even for Diophantine numbers 
with large partial quotients becomes substantially hard. So the question of establishing a Bryuno 
interpolation formula for e c (co), and obtaining the correct critical exponent fi if such a formula is 
indeed established, is still quite open. 
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In this paper we use Greene's method to compute the critical function when the distance of the 
rotation numbers from the resonances is of order 10 5 . As the computations close to resonances 
become very time-consuming we look at only three resonances (0, 1/2 and 1/3). We then use the 
conjectured Holder-continuity property of the function log s c (co) + fiB(co) to derive the corrections to 
the asymptotic behavior of log e c (a>), so improving significantly the agreement of the data with the 
conjectured value of /? = 1. 

Of course the problem is not completely solved, even from the numerical point of view, for two 
reasons. The first is that we consider only three resonances, so that a more exhaustive investigation 
would be needed. The second is the aforementioned very special choice of the sequences of rotation 
numbers tending to the resonances that we have to use. Nevertheless we improve the results existing 
in literature by one order of magnitude both in the distance from the resonance and in the value of 6, 
finding further support for the conjectured Bryuno's interpolation formula for the critical function. 

2 . Greene ' s method 

The main tool we use to determine numerically the break-down thresholds for analytic invariant 
curves for the standard map is Greene's method, known also as residue criterion. We now recall the 
main properties of the periodic solutions of the standard map used to formulate Greene's method, and 
sketch briefly its foundations, referring to the original paper fTSll for more details. 

We also recall that in 1131 and [ 24 1 some theorems are proved that go some way in the direction of 
proving the validity of Greene's method, at least in special cases. While a full rigorous justification 
of its use has not yet been achieved, Greene's method is considered one of the most accurate way to 
compute the critical function £ c for the standard map. 

If co is a rational number, given as the irreducible fraction p/q, then Birkhoff theory 1101 applies; 
its consequences for maps like the standard map T £ are the following. If s - (unperturbed, linear 
case) then there are trivially invariant curves with rational rotation number p/q, such that every point 
on them is a fixed point of the iterated map T° £ q . As the perturbation is turned on, only 2kq, ieff, 
points survive as fixed points of the q-th iterate of the map T e . These correspond to an even number 
(2k) of periodic orbits of period q. Such orbits - that we call perturbative - are the ones which will 
be studied within a perturbative framework; a simple perturbative calculation (see e.g. 0) shows that 
for the standard map the even number of such periodic orbits is indeed just 2. 

Of course this does not mean that such orbits are the only periodic ones for the standard map, but 
they are those which are obtained by continuation (in e) from unperturbed ones. In other words such a 
scenario does not consider the new periodic orbits arising when the perturbation is switched on. If we 
pass to the plane K 2 and consider the map T* , then the situation can be clarified in the following way. 
When co is irrational and satisfies the Bryuno condition, then the invariant curve with rotation number 
co of the unperturbed map survives for small values of e, while an invariant curve with rational rotation 
number p/q is suddenly destroyed; instead, only two discrete invariant sets of points {(^p,J] ( p)}j € z, 
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I = 1,2 survive, such that 




?Z a =?r + 2xp, £=1,2- 



(2.1) 



By taking the quotient in the first variable by the group of discrete translations by multiples of 27r, we 
obviously get two periodic orbits of period q, on which the motion has rotation number p/q. 

In [ 7 1 it is also proved that, for small values of s, each such periodic orbit lies on an analytic curve 
- called a remnant of the rational invariant curve of the unperturbed map -, and for rational numbers 
which approximate a Bryuno number oj such remnants approximate the invariant curve with rotation 
number oj. 

The basic idea of Greene's method consists in relating the break-down of an invariant curve with 
the loss of stability of nearby perturbative periodic orbits. In practice, the hypothesis behind Greene's 
method is that, if e < e c (oj), then there is a sequence of stable perturbative periodic orbits with rotation 
numbers pk/qk', as e grows beyond s c {oj), these periodic orbits lose stability in the large k limit. 

The criterion can be formulated more precisely in the following way. Let {(xf\yf^)}1 =1 be a per- 
turbative periodic orbit with rotation number pk/qk, approximating the irrational rotation number oj. 
Let 7~k(s) be the trace of the tangent dynamics along the periodic orbit: 



Then the periodic orbit is stable if -2 < 7~k(s) < 2, unstable otherwise. For historical reasons, the 
criterion is usually formulated in term of the residue Kk{s) of the orbit, related to the above trace by 



Therefore in terms of the residue the orbit is stable if < ^(fi) < 1, unstable otherwise. We then 
track, for a fixed value of s, the residue of those perturbative periodic orbits with rotation numbers 
Pk/qk which are stable for e = 0; if the residue diverges as k — > oo, then s > e c (oj), while if the residue 
tends to then s < e c (oj). 

It is actually conjectured (see I2"3l0 that if s < s c (oj), then the residue ^(e) tends exponentially to 
zero as k — ^ oo, with a rate of decay proportional to the width of the analyticity strip of the conjugating 
function u(a, e, oj) on the complex a plane for the values of e and oj considered. So Greene's method 
can also be used also to provide numerical information on the analytic properties of u in a, assuming 
this conjecture. 

An interesting question is what actually happens to the residue at the critical function s c (a>). It was 
originally conjectured that for noble rotation numbers (that is, rotation numbers which are obtained by 
applying a modular transformation to the golden mean, so that their continued fraction expansion has 
a "tail" of l's) it tends to a limit value, which should be about 0.25. We present below some numerical 
results which show that generally the situation is more complicate, and that such limit residue ^oo(e) 
could be not only different for different classes of rotation numbers, but could also be non-existent, 
and relate the behavior of the sequence of residues Hjc(e) for a fixed value of s along the sequence 
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of perturbative periodic orbits of rotation numbers pulqk to the arithmetic properties of the rotation 
number. 

From the practical, computational standpoint, the implementation of Greene's method faces some 
challenges if we wish to use it near resonances. The first is that, if a> is near a resonance, then the 
qk become soon quite large, that is we have to find many long periodic orbits, which takes a lot of 
computer time. 

The second, hardest, challenge is more subtle. In fact, if it happens that the rotation number of a 
periodic orbit is p/q « p' /q', with g » q' (the typical situation arising when approximating irrational 
rotation numbers close to small-denominator rationals) then it appears numerically that the periodic 
orbit of rotation number p/q tends to consist in lots (q is supposed to be large) of points accumulating 
near the points making the periodic orbit of rotation number p' /q'. The consequences for the compu- 
tation of the residue are dire, as in this case the matrix in (12.2ft has two very large, opposite, nearly 
equal in absolute value diagonal elements, so that when computing the trace the real data cancels and 
one is left with just the numerical error. Note that using a low precision with Greene's method so 
close to resonances gives essentially noise instead of the residue, so we get no values at all for B. 
We choose a brute-force solution to this precision problem, which consists in increasing the number 
of digits in the calculations until some data is left when computing the trace. Empirically, this could 
mean that one has to use hundreds of digits of precision in computing (12.2ft numerically, therefore also 
the periodic orbits must be known with such a precision: considering that one easily needs periodic 
orbits of period in the range of several tens of thousands - we actually reach orbits of length of the 
order of 150000 -, the calculation of a single value of e c can require a great amount of computer time. 

3. Numerical results 

3.1. Rotation numbers close to 0. Consider rotation numbers a> n = l/(n + y) = [n, 1°°], with n e N: 
in table [2 we give the values of the Bryuno function and of the critical function for rotation numbers 
oj nk , with [n^] a finite increasing sequence. Note that we reach values of rotation number close more 
than 2 x 10 -5 to the resonance value (0 in this case), which corresponds to values of n up to 60000. 
By fitting y = — loge c (a»„ t ) as a linear function of x = B(a>, lk ), we obtain 

y = ax + b a -0.9705, b = -1.9553. (3.1) 

As we see the slope is close to (but different from) 1: the relative difference is about 3.0%. 

One also realizes that the slope of the line increases if we neglect the rotation numbers a>„ corre- 
sponding to smaller values of n: this suggests that, if we consider just pairs of successive rotation 
numbers and evaluate the slope of the line passing through them, then we obtain an increasing func- 
tion. This can be formulated more precisely as follows. For n, m € N define 

, loge c (a;„)-loge c (fa; m ) 

A(a) n ,0} m ) = — — - — — — , (3.2) 

B(oj„) - B(cj m ) 

which measures the slope a of the line 

- log E c ((di) = aB(a>) + b, (3.3) 
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passing through the points (B(aj n ), - log e c (tu„)) and (B(a> m ), - log £ c (a> m )). We set - A(a> nk+l , a> nk ). 
In tablefflwe give also the values of the slopes A#: as we noted A^ steadily increases. 

The value B = 1 is anyhow still far from being reached: at best, just considering the last value of A^ 
in table^we obtain a value whose relative difference from 1 is greater than 1%. Moreover, though the 
values of the slopes increase as n —> oo, the convergence to 1 is very slow. In the next section we shall 
provide a heuristic argument which allows to guess the correction to the asymptotic behavior and so 
try to extrapolate a better value of B; this applies also to the cases considered in the next subsections. 

In table|2we give the values of the Bryuno function and of the critical function for a finite sequence 
of rotation numbers io nk , with oj nk = l/in^ + 1/(20 + y)) = [n^, 20, 1°°]: such numbers tend to as the 
previously considered ones, and share with them, essentially, the same Diophantine properties, as they 
have the same "tail" of 1 's in their continued fraction expansion, with the only difference that there 
is a partial quotient 20 before such a "tail". The distance of the rotation numbers considered is up to 
2 x 10~ 4 from 0, i.e. an order less than in the previous case: this is due to the fact the partial quotients 
go faster, and it becomes longer for the residue to reach the asymptotic value (so that periodic orbits 
with larger periods should be considered in order to obtain for the rotation numbers the same distance 
from the resonance value). 

As one can see, the values of the Bryuno function and of the critical function are comparable with 
those listed in table [fl the introduction of a larger partial quotient does not introduce any relevant 
change. As a consequence, also the slopes At, defined as before with the new definition of a>„ k , are 
very similar (as a look at the last column of table El immediately confirms). 

Note however that to compute numerically the critical function for rotation numbers of the form 
[rik, 20, 1°°] for given k is much more time consuming, since, in general, to obtain a reliable precision 
we are forced to reach periodic orbits with very high periods (say more than a hundred thousand), 
which requires a precision of about 600 digits. 

3.2. Rotation numbers close to 1/2. In table[3]we consider a sequence of rotation numbers tending 
to 1/2 of the form a>„ = 1/(2 + l/(n + y)) - [2,n, 1°°]. The rotation numbers considered are up to 
10" 5 close to the resonance value 1/2 (which correspond to values of n up to 20000). 

The fit for y = - log s c (a) nk ) as a linear function of x = B(a>, 1k ) gives 

y = ax + b a -0.9641, b = -1.6203. (3.4) 

Again we see the slope is not 1, and the relative error is now about 3.6%. It is greater than in the 
previous case because we stopped to smaller values of n; in fact the values of the slopes listed in table 
0show that again the function A^, defined exactly as before with the new definition for the rotation 
numbers u)„ k , is increasing in k. The relative difference from 1 of the last value of A^ is about 1.7%. 

3.3. Rotation numbers close to 1/3. In table|4]we consider a sequence of rotation numbers tending 
to 1/3 of the form a>„ = 1/(3 + l/(n + y)) - [3,n, 1°°]. The rotation numbers considered are up to 
5 X 10~ 6 close to the resonance value 1/3 (which correspond to values of n up to 20000). 
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The fit for y = - log e c {oj nk ) as a linear function of x = B(a) nk ) gives 

y = ax + b a - 0.9637, b = -1.6526. (3.5) 

while in the last column of tabled] we list the slopes A*, again defined as before with the new definition 
for the rotation numbers a), u ; the relative difference of the slope with respect to 1 is more than 3.7%, 
while the relative difference from 1 of the last value of is about 1.7%. 

3.4. Behavior of the critical residues and other rotation numbers. The behavior of the residue for 
e exactly equal to the critical function e c (a>) when the rotation number of the approximating periodic 
orbits tends to a) has been considered since the very first papers on the subjects (for example, in fTSll 
itself). In particular, one considers the sequence of residues ^(e c (w)) when k — » oo; it appears that 
this sequence has a limit only when u> is a number of so called "constant type", i.e. when to can be 

written as [ai a^, d°°]. This limit moreover seems to depend only on the integer d, and not at all 

from the "head" of the continued fraction expansion [a\,... ,ajf]. Unfortunately, a sound numerical 
evidence can be obtained only for J - 1 and for short "heads" in the continued fraction expansion, 
otherwise the partial quotients gets soon large and it becomes difficult to compute the critical residues 
with the accuracy required: therefore we state this more as a somewhat numerically founded and 
reasonable conjecture than else. In table|5]we give some values of the critical residue for a few values 
of d. 

If a rational number is not of constant type, then a limit does not seem to be achieved for the 
sequence of critical residues. In fact, it seems to happen that if co is a quadratic irrational, so that the 
sequence of the partial quotients a\ is eventually periodic, the sequence of critical residues is itself 
eventually periodic with the same period. In tables |6j 00 we can see the sequence of critical residues 
for some quadratic irrationals with short periods (resp. 2, 2 and 3). If the rotation number is not a 
quadratic irrational, so the partial quotients are aperiodic, the sequence of critical residues does not 
seem to have any regularity (but see below for a numerical difficulty). 

So far only quadratic irrational <x> have been considered. This is of course a limitation, due mainly 
to practical reasons; in fact, quadratic irrationals are the only irrationals with an eventually periodic 
continued fraction expansion, so they are particularly suited to Greene's residue criterion for two 
reasons: (1) the partial quotients are bounded, since they are periodic, so the approximants pklqk 
have denominators which do not grow too much and (2) if the period is reasonably small, one can 
tell whether the critical function has been reasonably approximated by looking at the sequence of 
residues over a span of periods and easily see whether it decreases or increases instead of being 
periodic. Instead, if the sequence of the partial quotients is aperiodic (and worst yet, unbounded) 
one can never be sure that the critical function has been obtained since the next periodic orbit to be 
considered (corresponding to the next approximant pklqk) could come from an abnormally high (or 
low) partial quotient a^. Note that in fl2l . where general irrationals are also considered in Subsection 
3.3, in the numerical calculations of the critical function, only the first ten partial quotients of the 
rotation numbers are retained, and all the others are set to 1, so one practically comes back to the case 
of noble numbers like ours. 
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This raises the question whether the algebraic (rather than just number-theoretic) properties of the 
rotation number have any role in the properties of the corresponding invariant curve. Lindstedt series 
expansion methods for instance do not care about the algebraic properties of oj, as the only relevant 
property is whether oj is a Bryuno number or not. "Phase space" renormalization group methods 
instead seem to work (or at least they have been applied) only in the case of quadratic irrational 
rotation numbers, so their results could depend on the algebraic layer. We expect that the algebraic 
properties of oj could show up, maybe, in discussing the smoothness of the natural boundary in the 
complex £ plane, but of course this is just a speculation. 

4. Discussion 

Despite intensive numerical calculations, the problem of confirming the conjecture expressed in 
dl . 10l > and estimating the exponent B cannot be considered completely settled even from the numerical 
point of view. In fact, as we noted earlier, only the three resonances 0/1, 1/2 and 1/3 have been 
considered, and only very special sequences of rotation numbers tending to such rational numbers 
have been used: considering other sequences of rotation numbers, in fact, means using numbers 
which have quite soon very large partial quotients, so that they are very bad from the numerical point 
of view. 

Moreover, a simple linear fit of log e c (oj) against B(a>), that is a fit which takes into account only 
the leading conjectured asymptotic behavior without any corrections, still gives results which are quite 
unsatisfying, as the difference between the estimated value of B and the conjectured value B = 1 is still 
of the order of a few percent. What is worst, the "running slopes" Ak defined in the preceding section 
continue to grow monotonically from below, slowly but steadily, so that one cannot even conclude 
that the conjecture is false or that the value of B is actually smaller than 1. Clearly corrections must 
be taken into account, or otherwise rotation numbers even closer to the resonances (and significantly 
such) must be considered, which is numerically unfeasible with current resources. 

Note also the apparently quite singular fact that for p(oj) the value 2 of the corresponding critical 
exponent seems to be obtained within a few percent much earlier. For instance for the rotation num- 
bers co close to 1 /2 listed in table [5J by using for the corresponding radii of convergence the values 
Pp{oj) computed by Pade approximants, we obtain for the slopes A' k = A'(oj nk+l , oj nk ), with 

v logp(^„)-logp(^,„) 

A ((On, 0J m ) = — - — - — , (4. 1) 

B((o n ) - B(oj m ) 

the values in the last column of table [5J Analogously for the rotation numbers oj close to 1/3 listed in 
tableEB again by using the values pp(co) computed by Pade approximants for the corresponding radii 
of convergence, we obtain the slopes A', in the last column of the same table. 

In ngureQ]we represent the analyticity domains for oj = [3,20, 1°°], [3,50, 1°°], [3, 100, 1°°] and 
[3,200, 1°°] as given by the poles of the Pade approximants [240/240]. As noted in Section ^ for 
oj getting closer to 1 /2 the poles tend to accumulate near the strongest singularity: therefore Pade 
approximants are not suitable for determining the critical function, but they can be fruitfully used in 
order to detect the radius of convergence. 
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Figure 1: Poles of the Pade approximant [240/240] for a> = [3, n, 1°°] and a = 1. 



The relative errors with respect to 2 for the values corresponding ton = 40, taken from tables[5]and 
ITOl are about 1.4% and 3.0%, respectively, therefore they are comparable with the errors for the last 
entries of the corresponding tables and 0] for the critical function for rotation numbers much closer 
to the resonance values: in the latter case indeed such errors are about 1.7%. And for larger values of 
n the relative errors become much smaller: for instance, for n = 100 and n = 200, we find from table 
[T0l errors about 1.0% and 0.5%, respectively. 

Of course it would be also interesting to have the slopes for the rotations numbers appearing in 
tables |3] and |4] To obtain the values of p(a>) numerically can be as hard as to determine the critical 
function e c (ll>). Also using the method of Pade approximants can be delicate, as in order to obtain 
reliable results a very high precision could be necessary. One could think of using the complex 
extension of Greene's method envisaged in flol . and the analysis, at best, could be as delicate as 
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in the present paper, where real values of e have been studied. We have also two more difficulties 
with respect to the case of s c (oj)- First one has to guess the direction in the complex plane where 
the singularities of the boundary of the analyticity domain are the closest to the origin; in this respect 
the results of [2| suggest, as a natural Ansatz, that, for rotation numbers close to p/q, they such 
singularities lie along the directions of the 2qth roots of - 1 . Next, for fixed a> close to a resonance 
value, the value of p(oS) should be much smaller than the value of s c (a>), again as a byproduct of the 
numerical analysis of [2] (and also that of (TTi). so that the value of p(a>) is expected to be harder to 
detect than e c (a>), as it should require more precision and hence more computing time. However we 
prefer to avoid any technical difficulties and to circumvent the problem by using the heuristic formula 
introduced in |;2j, say p(co) ~ pi(to), with 



where 77 = \a> - p/q\ if a> is close to the resonance p/q, C p / q is the numerical constant introduced in 
(one has C /i - 1, C 1/2 = -1/8 and C1/3 - -1/24), and Ac = 4n 2 x 0.827524 « 32.669338. 

In El we have already seen that there is a good agreement between the value pp(a>) of the radius 
of convergence found by Pade approximants and the value pi(a>) predicted by the formula (14.2ft . 
Furthermore the formula (14.2b becomes more and more reliable as a> approaches an rational value. 
See for instance tables |9] and which show how the difference between the two values pp(co) and 
p\{co) tend to shrink to zero when making the rotation number a> closer to the rational values 1/2 
and 1/3, respectively. So we can expect that the approximation we make by evaluating the radius of 
convergence p(co) with pi(cu) is very good for values much closer to the resonance values, as the ones 
we have considered are. 

Then we obtain the values listed in tables ^2 ^] and El for values of a> close, respectively, to 0, 
1/2 and 1/3 (the same for which we determined numerically the critical function); the slopes A' k are 
listed in the last columns of these tables. Of course, if we use the formula (I4.2I) . a slope approximately 
equal to 2 is expected, by the definition itself of p\{oS). The important fact is, in any case, that the 
discrepancy with respect to the value B = 2 (which in such a case is known to be the right one) is much 
smaller. In other words the asymptotic formula dl.9t is reached much earlier than the one which is 
believed to hold for the critical function. 

This different speed in reaching the asymptotic behavior of p{oS) and s c (a>) can be explained in 
terms of different corrections to the leading order when co —> p/q (and therefore B(a>) —> 00). We 
shall now try to compute such correction, at least heuristically, both for p{oS) and for e c (o>), and try to 
use them to extrapolate a better value of B. 

For what concerns p(a>), we shall assume the validity of the heuristic formula (I4.2t : this of course 
can introduce further corrections not accounted for, which we neglect assuming that they are smaller. 

Consider for example a>„ = l/(n + y), n —> 00. Then q = to n in (14.2ft and 




(4.2) 



B(o) n ) - - log o)„ + co n B(\/co n ), 
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which implies that log(/i + y) - B(a> n ) - B(y)/(n + y). Therefore in first approximation we have that 

log(« + 7) a B(a> n ) - B(y)e- B(0J "\ 

as the leading behavior of B(oj n ) is just log(n + y). This gives 

log Pl (w n ) * logdCo/il" 1 ^) - 2fi(w„) + 2B(y)e' BM , (4.3) 

that is the correction to the linearly growing asymptotic behavior is exponentially small. An analogous, 
slightly more complicated, computation for a> n - 1 /(q + 1 /(n + y)) gives a coiTection of the form 
B(u> n ) exp(-qB(tL> n )), that is still essentially exponentially small. This explains the exceptional rapidity 
of the approach to the scaling behavior for p(a>). 

Quantitatively, a fit of the numerical data of table [fusing <I4.3I > to model the data gives for B the 
value 1.9999989, whose difference from the correct value of 2 is of the order of 1(T 7 , while a straight 
linear fit gives 2.00091, whose error is three orders of magnitude larger. 

Analogously, a fit of the numerical data of table [T3*1 using 

log Pl (^„) * logOCc/iP 1 ^) -8B(a) n ) + {b + cBico))- 3 ^ , (4.4) 

gives for B the value 2.000000287, whose difference from the correct value of 2 is of the order of 10~ 7 , 
while a straight linear fit gives 1 .99984, whose error is again three orders of magnitude larger. Also a 
comparison between the mean-square distances of the data from the corresponding fits is remarkable: 
we obtain 2.389 x 10~ 8 for the fit by (§0), and 0.0000858 for the linear fit. Note that also fits with 
either b = or c - in (14.4ft are worse: for b - we obtain a value B - 1.99966 (with mean-square 
distance 0.0000413), while for c = we obtain a value B = 1.99967 with mean-square distance 
0.0000393). 

To compute the correction to the leading behavior of s c (o») is of course quite another matter, since 
we don't even have a proof or at least a very strong theoretical argument for the leading order. So the 
following argument is more a qualitative explanation rather than a quantitative attempt to extrapolate 
seriously the value of B (but we shall try nevertheless). 

As before we shall consider only the case a> n = l/(q + l/(n + y)), n — > oo, and we shall set 
q n = \u) n - llq\. Let 

log E c (o) n ) + B(aj n ) = C(oj n ), (4.5) 

where the function C(a») is believed to be continuous (see [26 28 1). Let then c = lim w _>o C(oj). It is 
also conjectured (see E71 ) that C(o>) is Holder-continuous with some exponent a (in the quoted paper 
it is suggested that a could be 1/2): so, by recalling that n„ « e -i B ^n) q, a reasonable guess in 
\4.5\ could be 

log e c (co n ) - (const.) - BB(co n ) + 0{e- aqBM ). (4.6) 
Despite the rough, qualitative nature of the argument above, we can try to fit the data with the formula 
<I4.6I > and see whether the growth of the slopes is such that the value of 1 can actually be reached. 
The "best" value of a is obtained by choosing it in such a way that the mean square distance of the 
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experimental data from the values obtained from the fit is minimal. As an alternative, we performed 
also nonlinear fits using Levenberg-Marquardt method (see I2"T1|2"9"1). obtaining consistent results. 

Fitting the data relative to £ Q (co n ) for the sequence considered in table[2with the formula (14.6b . we 
obtain 

loge c (w„) ~ -2.34630 + 1.00359 B(u n ) + 1 .59684 e" 033025 ^, (4.7) 
that is finally a value much closer to 1 than the straight linear fit, which gave 0.97052. Moreover, the 
mean-square distance of the data from the fit is 0.000210 in the case of the fit with corrections, while 
is much larger, that is 0.0396, in the case of the linear fit (see figure Et). 




3.2 3.4 3.6 3.8 4 4.2 4.4 



(c) Sequence co„ k listed in table|4] 

Figure 2: Numerical values of - log£ c (<y„ ( ), obtained with Greene's method, versus B(a>„ k ) for the 
sequence oj„ k listed in tablesQ[3] and[3J the error bars are less than the size of the points. The solid 
curve corresponds to the fits 14.71 . 14.81 and 14.91 . 

If we consider the case a> — > 1/2 as in subsection !3.2l the fit with formula (I4.6I > gives 

loge c (w„) a -1.86364 + 1.00308 B(a) n ) + 1 .43766 e -069671 ^ (4.8) 

with mean-square distance d = 0.0000512 (the linear fit would give ft = 0.96413 and d = 0.0124); in 
figure EJ) we plot the data together with the fit. 
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Finally, in the case a> —> 1/3 as in subsection !3.3l the fit with formula ( 14.61 ) gives 

log£ c (w„) * -1.84393 + 1. 00344 B(a) n ) + 1.82643 e~ 1 0300 B{0J "\ (4.9) 

with mean-square distance 0.0000403 (the linear fit would give B = 0.96369 and d = 0.00832); in 
figure 0; we plot the data together with the fit. 

The results, together with the interpolation formula (14.61) . hint at a value of a close to 1/3, while 
C p (oj), according to the formula (14 .4t . seems to be Holder-continuous with any exponent a p < 1 (and 
a p = 1 in 0). So while in both cases the corrections are exponentially small in B(a>), the coefficient 
in the exponential is about three times larger for log p(a»), leading to smaller corrections and faster 
approach to the asymptotic regime. 

If we try to plot C p (cS) for the values of a> close to listed in tabled and use the values of p(a>) in 
table ^2 we And the behavior represented in figure^, which also support the smoothness conjecture. 
Analogously, if we plot C(co) for the same set of values of a>, by using the values of e(oj) listed in 
table [0 we find the behavior in figure 0). 

In figures |3j; and |3ji we represent the functions C p (oj) and C(a>) for the values of od close to 1 /2 
listed in table |3j and in figures |3^ and|3f we represent the functions C p {u) and C(oj) for the values of 
a) close to 1 /3 listed in table @] 

While the variation in the case of £ c (u>) is larger, all plots support the conjecture of a function 
which is not only bounded but also Holder continuous close to the resonances, but of course a deeper 
numerical investigation is needed in order to draw more quantitative deductions. 

To conclude, we note also that the value 6 = 2, which holds for the radius of convergence, is found 
both in fTTl and in flD for the critical function, if the computations are made without requiring a 
high precision (this is attributed to the low precision in ifTTl and to the truncations and approximations 
due to the numerical implementation of the renormalization group method in [12]). We find this 
phenomenon at least very curious: it would be in fact quite interesting to understand why truncation 
and approximation errors in the numerical computations give a different value of B (and exactly the 
one holding for the radius of convergence) instead than just gibberish. 

5. Conclusions 

We conclude by some general remarks about the advances which have been made and the conclu- 
sions which can be drawn from our analysis. 

(1) The numerical results of 1111 have been improved by an order of magnitude, both in the size 
of Aco and in the order of the numerical errors. Moreover, an heuristic argument providing 
corrections to the leading order has been given: the analysis of the numerical data, taking 
into account the conjectured form of the corrections, supports both (11.101) with B = 1 and the 
continuity of the function C(a>). A stronger support would require getting closer to the reso- 
nances, and considering more resonances rather than just 0/1, 1/2 and 1/3: all these actions 
are clearly feasible but would require significantly more computer time, which has already 
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Figure 3: Plot of C p (co) = logp(co) + 2B(a>) and C(co) = log£ c (<jj) + B(co) versus to for the sequence 
co„ t listed in tablesQ(plots a and b),[3](plots c and d) and0(plots e and f). 



reached the order of several CPU years on Compaq Alpha computers for the calculations of 
the present paper. 

(2) While a reasonably complete analysis of what happens for sequences of rotation numbers 
which are not the best ones cannot be practically done, the study of a sequence like [n, 20, 1°°] 
shows that the behavior of the critical function along this sequence is the same that along the 
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"best" sequence [n, 1°°]. It is likely that this holds, in the limit of long periodic orbits, at least 
for all sequences of noble numbers tending to a rational value. We could investigate very 
few non-noble sequences and no sequence at all made by something different than quadratic 
irrationals (which have measure 0). For a sequence like [n, 2°°] the values of e c (a>) seem to be 
comparable to the ones of the "best" sequence quoted above. 

(3) The study of the functions C(u>) = log s c (a>) + B(a>) and C p (a>) = logp(w) + 2B(a>) seem to 
suggest that they depend smoothly on a>. In general such functions look as continuous in their 
arguments, as also the comparison between the two sequences [n, I 00 ] and [20, n, l°°] seem 
to support. This is in contrast with the conclusions made in l22ll . where doubts were raised 
about continuity of the function log£ c (tL>) + BB{oS) for any choice of 6. 

(4) Some interesting conclusions can be drawn for the behavior of the critical residues. It appears 
that a limit value is obtained only when a> is of constant type, and this limit seems to depend 
only on the "tail" of the expansion. If od is not of that form but still a quadratic irrational, 
then the sequence of the partial quotients is eventually periodic. In this case it appears that 
"Rkisdoj)), for large k, approaches a periodic sequence of values with the same period of the 
partial quotients. If a> is not a quadratic irrational, then it is difficult to draw any conclusion 
at all. If the partial quotients are bounded, then the critical residues seem to be bounded away 
from zero and oscillating in an apparent random way, but bounded; for unbounded sequences 
of partial quotients no numerical data at all could be obtained. We note that such a scenario 
is consistent with that arising within the renormalization group approach as described, for 
instance, in llT2l . 
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Tables 

Table 1: Values of the Bryuno function, of the critical function and of the running slopes A k = 
A{u) nt ,u)„ k j) corresponding to a finite sequence of rotation numbers a)„ k = l/(n k + y) = [n k , 1°°]. 
The error on e c (u>„ t ) is of 1 unit on the last digit, and the corresponding slopes are computed with 
consistent accuracy. 
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6 


[7000, r°] 


8.85393 


0.0013265 


0.9724±0.0002 


7 


[10000, 1°°] 


9.21053 


0.00093627 


0.9770±0.0002 


8 


[12000, 1°°] 


9.39284 


0.00078320 


0.9793±0.0001 


9 


[15000, 1°°] 


9.61593 


0.00062927 


0.9808±0.0001 


10 


[18000, 1°°] 


9.79823 


0.00052610 


0.9823±0.0002 


11 


[20000, 1°°] 


9.90358 


0.00047433 


0.9833+0.0004 


12 


[25000, 1°°] 


10.12671 


0.00038081 


0.9842+0.0002 


13 


[30000, 1°°] 


10.30902 


0.00031816 


0.9859+0.0003 


14 


[40000, 1°°] 


10.59668 


0.00023955 


0.9865+0.0003 


15 


[50000, 1°°] 


10.81982 


0.000192161 


0.9879+0.0003 


16 


[60000, 1°°] 


11.00213 


0.000160443 


0.9895+0.0002 



Table 2: Values of the Bryuno function, of the critical function and of the running slopes A k = 
A(a) nk ,LJ„ k j) corresponding to a finite sequence of rotation numbers a>„ k = \j(n k + 1/(20 + y)) = 
[n k ,20, 1°°]. The error on e c (a>„ k ) is of 1 unit on the last digit, and the corresponding slopes are 
computed with consistent accuracy. 



k 




B((o nk ) 


£ c (oJ,i k ) 


A k 


1 


[500, 20, 1°°] 


6.22088 


0.016303 




2 


[700,20, 1°°] 


6.55556 


0.011926 


0.9341+0.0004 


3 


[1000,20, 1°°] 


6.91089 


0.008535 


0.9415+0.0006 


4 


[2000,20, 1 00 ] 


7.60247 


0.004421 


0.9512+0.0005 


5 


[4000,20,1 00 ] 


8.29483 


0.002271 


0.962 ±0.001 
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Table 3: Values of the Bryuno function, of the critical function and of the running slopes A k = 
A{co nt ,(i) nk j) corresponding to a finite sequence of rotation numbers oj, lk = 1/(2 + l/(n k + y)) = 
[2, n k , 1™]. The error on s c (cj nk ) is of 1 unit on the last digit, and the corresponding quantities are 
computed with consistent accuracy. 



k 


CO,,, 


B(a> n . ) 


£c(^m ) 


A k 


1 


[2,500, 1°°] 


3.80022 


0.12872 




2 


[2,700, 1°°] 


3.96840 


0.109967 


0.9362±0.0005 


3 


[2, 1000, 1°°] 


4.14674 


0.092932 


0.9438±0.0001 


4 


[2,2000, 1°°] 


4.49337 


0.066777 


0.9535±0.0001 


5 


[2,4000, 1°°] 


4.84001 


0.047805 


0.9642±0.0001 


6 


[2,7000, 1 00 ] 


5.11987 


0.036420 


0.9720±0.0002 


7 


[2, 10000, 1°°] 


5.29823 


0.030598 


0.9766±0.0003 


8 


[2, 13000, 1°°] 


5.42943 


0.026909 


0.9792±0.0006 


9 


[2, 17000, 1°°] 


5.56357 


0.023591 


0.9810±0.0006 


10 


[2,20000, 1 00 ] 


5.64484 


0.021780 


0.983 ±0.001 



Table 4: Values of the Bryuno function, of the critical function and of the running slopes A k = 
A{u) nk ,(i) nt _,) corresponding to a finite sequence of rotation numbers io„ t = 1/(3 + l/Oi* + y)) = 
[3, n k , 1™]. The error on s c (cj nk ) is of 1 unit on the last digit, and the corresponding quantities are 
computed with consistent accuracy. 



k 




B(a>n k ) 




A k 


1 


[3,500, 1°°] 


3.17069 


0.244787 




2 


[3,700, 1°°] 


3.28264 


0.22044 


0.9358±0.0001 


3 


[3, 1000, 1°°] 


3.40139 


0.197080 


0.9433±0.0001 


4 


[3,2000, 1°°] 


3.63230 


0.158153 


0.9529±0.0001 


5 


[3,4000, 1°°] 


3.86330 


0.126588 


0.9637±0.0001 


6 


[3,7000, 1°°] 


4.04983 


0.105608 


0.9715±0.0001 


7 


[3, 10000, 1°°] 


4.16872 


0.094035 


0.9763±0.0002 


8 


[3, 13000, 1°°] 


4.25617 


0.086319 


0.9787±0.0003 


9 


[3, 17000, 1°°] 


4.34559 


0.079072 


0.9807±0.0003 


10 


[3,20000, 1°°] 


4.39977 


0.074973 


0.9826±0.0005 
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Table 5: Critical residues for some rotation numbers oj. The error on s c (a)„ k ) is of 1 unit 

the last digit. 





s c (u>) 




[n 


0.971635406 


0.250088 


[2°°] 


0.957445408 


0.2275138 


[3°°] 


0.890863502 


0.202230 


[4°°] 


0.80472544 


0.17923 


[10, 2°°] 


0.481985986 


0.22751 


[1,3,2°°] 


0.829500533 


0.22751 


[7,3°°] 


0.615071885 


0.2022 


[1,2,4°°] 


0.86423037 


0.1792 



Table 6: Residues of critical periodic orbits for to = V3 - 1 = [1, 2, 1, 2, 1, 2, 



s c (a>) 


0.876067426 


approximant 


residue 


3/4 


0.24871 


8/11 


0.18612 


11/15 


0.25216 


30/41 


0.18516 


41/56 


0.25275 


112/153 


0.18493 


153/209 


0.25288 


418/571 


0.18487 





0.876067426 


approximant 


residue 


571/780 


0.25291 


1560/2131 


0.18486 


2131/2911 


0.25292 


5822/7953 


0.18485 


7953/10864 


0.25292 


21728/29681 


0.18485 


29681/40545 


0.25292 


81090/110771 


0.18486 
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Table 7: Residues of critical periodic orbits for a> = ( V3 - l)/2 = [2, 1, 2, 1, 2, 1, 



e c {oj) 0.9402827 

approximant residue 

3/8 0.19574 

4/11 0.24746 

11/30 0.18763 

15/41 0.25145 

41/112 0.18556 

56/153 0.25254 

153/418 0.18503 

209/571 0.25282 



e c (w) 0.9402827 

approximant residue 

571/1560 0.18490 

780/2131 0.25290 

2131/5822 0.18486 

2911/7953 0.25292 

7953/21728 0.18486 

10864/29681 0.25292 

29681/81090 0.18486 

40545/110771 0.25293 



Table 8: Residues of critical periodic 



e c (oj) 


0.9402827 


approximant 


residue 


3/5 


0.2242 


4/7 


0.2639 


7/12 


0.2278 


18/31 


0.2222 


25/43 


0.2660 


43/74 


0.2270 


111/191 


0.2227 


154/265 


0.2656 


265/456 


0.2272 



forw= V572- 1 = [1,1,2,1,1,2,1,1,2. 



s c {oj) 


0.9402827 


approximant 


residue 


684/1177 


0.2226 


949/1633 


0.2656 


1633/2810 


0.2271 


4215/7253 


0.2227 


5848/10063 


0.2656 


10063/17316 


0.2271 


25974/44695 


0.2227 


36037/62011 


0.2656 


6201 1/106706 


0.2272 
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Table 9: Radius of convergence for some values of the rotation number to close to 1/2 and slopes 
A' k = A'(lo n , , w„ t ] ). The value pi(toi) is given by the formula <4.2> . while pp(oj) is the value obtained 
numerically by using Pade approximants. the two values for the slopes correspond to the values 
pi(a>) and pp(oS), respectively. One has rj = \a> - 1/2]. 



k 


(o„ k 


*7 


Pi (AW.) 


Pp(oj nk ) 


K 


1 


[2, 10, l 00 ] 


0.0224860 


0.51409 


0.51052 




2 


[2, 12, 1°°] 


0.0190577 


0.43571 


0.43355 


2.19667/2.17013 


3 


[2, 15, 1°°] 


0.0155106 


0.35462 


0.35352 


2.14426/2.12464 


4 


[2,20,1°°] 


0.0118382 


0.27066 


0.27024 


2.09658/2.08449 


5 


[2,30,1°°] 


0.0080339 


0.18368 


0.18361 


2.05439/2.04822 


6 


[2,40,1°°] 


0.0060801 


0.13901 


0.13902 


2.02821/2.02484 


6 


[2,50,1°°] 


0.0048906 


0.11181 


0.11184 


2.01612/2.01480 



Table 10: Radius of convergence for some values of the rotation number to close to 1/3 and slopes 
A', = A'(a)„ k , w„,_j ). The value pi(a>) is given by the formula <4.2t . while p/>(w) is the value obtained 
numerically by using Pade approximants; the two values for the slopes correspond to the values 
pi(ai) and pp(cj), respectively. One has rj = \a> — 1/3]. 



k 






Pi(oj nk ) 




K 


1 


[3, 10, 1°°] 


0.0101459 


0.62329 


0.61993 




2 


[3, 12, 1°°] 


0.0085791 


0.55734 


0.55524 


2.28295/2.24934 


3 


[3, 13, 1°°] 


0.0079642 


0.53038 


0.52858 


2.23762/2.22067 


4 


[3,20,1°°] 


0.0053033 


0.40444 


0.40400 


2.17212/2.15360 


5 


[3,30,1°°] 


0.0035899 


0.31180 


0.31182 


2.09982/2.09051 


6 


[3,40,1°°] 


0.0027132 


0.25871 


0.25872 


2.06311/2.06339 


7 


[3,50,1°°] 


0.0021807 


0.22364 


0.22360 


2.04490/2.04795 


8 


[3, 100, 1°°] 


0.0011006 


0.14177 


0.14179 


2.02455/2.02313 


9 


[3,200,1°°] 


0.0005529 


0.08959 


0.08961 


2.00902/2.00866 
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Table 11: Values of the radius of convergence and of the slopes A' k = A'(p)„ k ,co nk ) corresponding 
to a finite sequence of rotation numbers a>„ k = l/(n k + y) = [n k , 1°°]. The radius of convergence is 
computed with the formula <4.2t . 



K 






A' 

A k 


1 


[500, 1°°] 


0.000130355 




2 


[700, 1°°] 


0.0000665545 


2.0042837 


3 


[1000, 1°°] 


0.000032629 


2.0030298 


4 


[2000, 1°°] 


0.00000816229 


2.0018183 


5 


[4000, 1°°] 


0.0000020412 


2.0009090 


6 


[7000, 1°°] 


0.000000666603 


2.0004825 


7 


[10000, 1°°] 


0.000000326653 


2.0003028 


8 


[12000, 1°°] 


0.000000226847 


2.0002303 


9 


[15000, 1°°] 


0.000000145185 


2.0001882 


10 


[18000, 1°°] 


0.000000100824 


2.0001536 


11 


[20000, 1°°] 


0.0000000816683 


2.0001329 


12 


[25000, 1°°] 


0.0000000522684 


2.0001129 


13 


[30000, 1°°] 


0.0000000362978 


2.0000921 


14 


[40000, 1°°] 


0.0000000204177 


2.0000730 


15 


[50000, 1°°] 


0.0000000130674 


2.0000565 



Table 12: Values of the radius of convergence and of the slopes A' k = A'(a>„ k ,a)„ tl ) corresponding 
to a finite sequence of rotation numbers co„ k = 1/(2 + l/(n k + y)) = [2,n k , 1°°]. The radius of 
convergence is computed with the formula 14.21 . 



k 


U>n k 




K 


1 


[2,500,1°°] 


0.011405915 




2 


[2,700, 1°°] 


0.008152279 


1.9968638 


3 


[2, 1000, 1°°] 


0.005709327 


1.9973651 


4 


[2,2000,1°°] 


0.002856258 


1.9980597 


5 


[2,4000, 1°°] 


0.001428528 


1.9987793 


6 


[2,7000,1°°] 


0.000816400 


1.9992292 


7 


[2, 10000, 1°°] 


0.000571507 


1.9994593 


8 


[2, 13000, 1°°] 


0.000439632 


1.9995765 


9 


[2, 17000, 1°°] 


0.000336196 


1.9996572 


10 


[2,20000,1°°] 


0.000285770 


1.9997123 
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Table 13: Values of the radius of convergence and of the slopes A' k = A'(a)„ k ,a)„ k ) corresponding 
to a finite sequence of rotation numbers u n . = 1 /(3 + 1 /(n k + y)) = [3, n k , 1°°]. The error on e c (a) nk ) 
is of 1 unit on the last digit, and the corresponding logarithm is computed with consistent accuracy. 



/V 


0)„ k 




A' 

A k 


1 


[3,500, 1°°] 


0.04873028 




2 


[3,700, 1°°] 


0.03895268 


2.0004598 


3 


[3, 1000, 1°°] 


0.03071760 


2.0000489 


4 


[3,2000, 1 00 ] 


0.01935701 


1.9997910 


5 


[3,4000, 1°°] 


0.01219611 


1.9997289 


6 


[3,7000, 1°°] 


0.00839894 


1.99977 '44 


7 


[3, 10000, 1°°] 


0.00662168 


1.9998204 


8 


[3, 13000, 1°°] 


0.00555920 


1.9998501 


9 


[3, 17000, 1°°] 


0.00464887 


1.9998731 


10 


[3,20000, 1°°] 


0.00417153 


1.9998900 
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